Detailed /3-decay study of 33 Ar 
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The proton-rich nucleus 33 Ar has been studied by detailed proton and 7-ray spectroscopy at the 
low-energy facility of SPIRAL at GANIL. Proton and 7-ray singles and coincidence measurements 
allowed to establish a quasi complete decay scheme of this nucleus. By comparing the proton 
intensity to different daughter states, tentative spin assignments have been made for some of the 
states of 33 CI. The Gamow- Teller strength distribution is deduced and compared to shell-model 
calculations and a quenching factor is determined. States close to the isobaric analogue state are 
searched for with respect to isospin mixing. 



PACS numbers: 
spectroscopy 



21.10.-k Properties of nuclei, 23.50.+Z Decay by proton emission, 29.30.Ep Charged-particle 



I. INTRODUCTION 

An important part of information about the struc- 
ture of the atomic nucleus is obtained by nuclear f3 de- 
cay. Studies close to the valley of nuclear stability en- 
able to precisely test our understanding of this structure, 
whereas experiments and theoretical work with nuclei far- 
away from the stability line enables investigations of the 
evolution of nuclear structure as a function of isospin. 
However, further away from stability and in particular 
on the proton-rich side of the valley of stability, nuclei no 
longer decay by simple j3 decay, but by /3-delayed parti- 
cle emission [l| . Therefore, spectroscopic studies become 
more difficult, as 7-ray spectroscopy has to be combined 
with particle detection, i.e. proton and a-particle detec- 
tion. 

Gamma detection is commonly achieved with high- 
efhciency germanium detectors which allow to resolve dif- 
ferent 7 lines with good precision. Charged particles are 
detected with silicon detectors which can reach modest 
energy resolutions of about 10-20 keV. A combination of 
charged-particle detection and 7-ray measurements may 
permit to establish even complicated /3-decay schemes 
and compare them to theoretical predictions from micro- 
scopic models as the nuclear shell model. 

A long-standing problem is the observation of " quench- 
ing" of the Gamow- Teller (GT) strength, when exper- 
imental results, either from nuclear j3 decay or from 
charge-exchange reactions, are compared with theoretical 
results. Having a slight mass dependence this quenching 
reaches a value of about 0.6 for the transition strength. 
The quenching is seen from a comparison of the ex- 
perimental and the theoretically calculated transition 
strength B(GT) 0. 

Although this quenching is known since many years 



and has been observed for basically all nuclei where these 
kinds of studies has been performed, two possible expla- 
nations are still discussed in the literature (see 0): i) 
problems in the shell-model calculations which do not 
take into account enough orbitals to describe the tran- 
sition strength correctly. In particular, intruder orbitals 
giving rise to higher-energy excitations are usually not 
included to keep the model calculations feasible, ii) Sub- 
nuclear excitations, especially of the A resonance, may 
shift transition strength to excitation energies as high as 
100 MeV and are therefore rather difficult to observe ex- 
perimentally. A paper of Arima |3| retracing to some 
extent the history of this quenching shows that the GT 
strength quenching is to a large extent due to two-particle 
two-hole excitations not included in most calculations. 

Another question of debate is the isospin mixing or 
isospin purity of nuclear states. In the ft decay of proton- 
rich nuclei, the Fermi transition populates the isobaric 
analogue state (IAS) of the decaying ground state. For 
/3-emitting nuclei with isospin projection T z < 0, this 
IAS is within the Qec window and decays either by 7 
emission or, for sufficiently proton-rich nuclei, by proton 
emission. However, as this IAS has an isospin quantum 
number equal to the one of the parent nucleus, proton 
emission from this state to states in the /3-delayed proton 
(/?p) emission daughter nucleus is forbidden by isospin 
conservation and thus only permitted by a small isospin 
impurity of either the IAS or of the final state after pro- 
ton emission. As the IAS usually lies in a region with a 
much higher density of nuclear states than the daughter 
state, often a nuclear ground state, the IAS has gener- 
ally a higher degree of isospin impurity than its daughter 
state. If this is the case, the isospin mixing has its origin 
mainly in a mixing of the IAS with nearby lying states 
having a one-unit lower isospin quantum number, but 
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the same spin quantum number. According to pertur- 
bation theory, the mixing gets larger with decreasing en- 
ergy difference between the two states. Of course, mixing 
is possible not only with one state, but with any state, 
however, with rapidly decreasing amplitude as the energy 
difference increases. 

The decay of 33 Ar, the subject of the present work, 
was studied several times in the past (see P4TT|). Our 
experimental data are compared to the most recent work 
from these authors and experimental averages are con- 
fronted to theoretical shell-model calculations using dif- 
ferent effective interactions optimized for the 2sld nu- 
clear shell-model space. In particular, we will determine 
the Gamow- Teller strength distribution and investigate 
the possibility to identify the state(s) responsible for 
isospin mixing in 33 CI. The isospin mixing of the IAS 
was suggested to be as much as 20% for 33 CI 0], a huge 
value compared to "today" standards. 

Hardy et al. Q and Borge et al. [§[ studied the decay 
of 33 Ar in detail in order to extract the GT strength 
distribution. Schardt et al. [|| and Garcia et al. [ll| also 
determined limits of possible scalar contributions to (3 
decay and the Fermi-to-Gamow- Teller ratio of the decay 
of 33 Ar. We will in particular use this ratio in our analysis 
and comparison to theory. 

In the present work, the use of high-resolution sil- 
icon surface-barrier detectors and of a high-efficiency 
7-detection array permitted to determine the proton 
branches to the ground as well as to the first and second 
excited states in the proton daughter. In particular, the 
GT distribution in the full Qec window was established 
and a detailed comparison to shell-model calculations was 
performed. 



II. EXPERIMENTAL DETAILS 

The isotopes of interest were produced by pro- 
jectile fragmentation of a 36 Ar 18+ primary beam at 
95 MeV/nucleon with intensities between 4 /zA and 8 fiA 
provided by the coupled cyclotrons of GANIL. This beam 
was fragmented in the SPIRAL carbon target. The SPI- 
RAL ECR source then produced a low-energy beam of 
33 Ar 3+ , which was directed to the SPIRAL identifica- 
tion station, where the experimental setup was mounted. 
The beam line to this detection station was regularly op- 
timized by means of a stable 40 Ar 3+ beam produced with 
the ECR source. Secondary beam intensities of 1200 pps 
on average were obtained for 33 Ar 3+ . The 33 Ar beam was 
contaminated by small amounts of 35 Ar possibly trans- 
mitted as a gas to the detection setup. 

Beam profilers and Faraday cups close to the detection 
setup allowed to optimize the beam position in the exper- 
imental setup and to monitor the beam intensity (mainly 
stable isotopes or molecules). The secondary beam was 
finally intercepted in the center of the experimental setup 
by means of a 0.9 /xm thick aluminized Mylar foil (1 cm 
x 2 cm) mounted on a thin metallic frame. 




FIG. 1: (Color online) Schematic view of the Silicon Cube 
consisting of 6 DSSSDs backed by 6 PIPS detectors. The 
position of the mylar foil catcher and the entrance of the sec- 
ondary beam are also shown. 



The standard detection setup of the SPIRAL identi- 
fication station was replaced by the Silicon Cube detec- 
tor [r2j and three high-efficiency germanium detectors 
from the EXOGAM array [H|. The Silicon Cube con- 
sists of six double-sided silicon strip detectors (DSSSDs) 
with 16x16 strips and a pitch of 3 mm forming a cube 
surrounding the mylar catcher foil. These DSSSDs were 
backed by six large-area (50x50 mm 2 ) Passivated Im- 
planted Planar Silicon (PIPS) detectors used to detect 
f3 particles. The thicknesses of the different detectors are 
given in Table |TJ The arrangement of the different de- 
tectors can be seen on figure [1] A photo of the setup 
installed at the low-energy identification station of SPI- 
RAL at GANIL is shown in figured 



DSSSD - PIPS position 


1 


2 


3 


4 


5 


6 


DSSSD thickness (^m) 


300 


287 


270 


64 


1000 


288 


PIPS thickness (fim) 


300 


300 


500 


1473 


150 


1498 



TABLE I: Thicknesses of the different silicon detectors used 
in the present study. The arrangement of the detectors can 
be seen in figure [T] 

The 192 signal channels from the DSSSDs were read- 
out via a printed circuit board incorporated in the de- 
tector housing. The 16-channel pre-amplifier cards were 
mounted directly on the detector chamber enabling thus 
a rather compact configuration. The six channels of the 
PIPS detectors were connected to their pre-amplifiers via 
standard LEMO connectors. This setup reaches a geo- 
metrical efficiency for proton detection of 54% [l2[ • 

The Silicon Cube was surrounded by three EXOGAM 
clover detectors. These detectors were calibrated in ef- 
ficiency and energy by means of standard calibration 
sources and yielded a total efficiency of about 3.2% at 
1332 keV. 

Due to technical problems, DSSSD 2 did not function 
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FIG. 2: (Color online) Photo of the setup as installed at the 
low-energy identification station of SPIRAL at GANIL. The 
silicon cube detector is seen in the center of the photo. On 
the left, the vacuum pipe to insert the catcher is visible. Two 
out of the three EXOGAM clover detectors are mounted and 
displayed on the bottom and right-hand side of the photo. 



correctly and was excluded from the present analysis. 
The thickness of detector 4 limited its use to low-energy 
protons. Therefore, it was also excluded from the anal- 
ysis. Due to its large thickness, DSSSD 5 yields a quite 
high /3 tail thus worsening its energy resolution. We used 
this detector only for the high-energy part of the spec- 
trum (above 6 MeV) where the other detectors did no 
longer fully stop the protons. Therefore, we used detec- 
tors 1,3, and 6 to determine proton energies and branch- 
ing ratios for proton lines below 6 MeV and detector 5 
above. The spectrum above 6 MeV was corrected for the 
geometrical efficiency of detector 5 with respect to de- 
tectors 1, 3, and 6. The /3 signal in the PIPS detectors 
was used in the analysis to reduce the /3-p summing by 
requiring to have the f3 signal in any of the other back 
detectors of the cube. This condition permitted to sig- 
nificantly reduce the /3 tails in the proton spectra. 

All silicon detectors were first calibrated by means of a 
triple a source ( 239 Pu, 241 Am, 244 Cm). The final energy 
calibration of the DSSSDs was performed using the /3- 
delayed proton peaks to the ground state of 32 S from the 
IAS recently remeasured [3] to be at 5547.8(9) keV and 
the 3/2+ state at 3971.9(12) keV [H excitation energy 
in 33 CI, respectively. An energy resolution of 50 keV was 
obtained for the DSSSD sum spectrum from detectors 1, 
3, and 6. Events in the DSSSDs were accepted only if the 
signals in the front and back side were within ±100 keV 
of each other. 

Different runs optimized for 33 Ar were performed yield- 
ing a total of 28 h for this isotope. Several runs were also 
performed for 31 Ar. These latter settings were strongly 
contaminated by 33 Ar and yielded rather low production 
rates for 31 Ar. These data were thus only used to study 
the performance of the Silicon Cube [l2[ . 



For the data taking with the different argon isotopes, 
only the DSSSDs were allowed to trigger the data ac- 
quisition. After a trigger, all channels were readout by 
means of the GANIL data acquisition system via VXI 
and VME modules. 



III. ANALYSIS PROCEDURE AND 
EXPERIMENTAL RESULTS 

The experimental data were analyzed as singles and 
coincidence data. The different peaks were fitted and in- 
tegrated by means of the PAW package from the CERN 
library. Most of these peaks are due to /3-delayed pro- 
ton emission from states in the /3-decay daughter nucleus 
33 C1 to the ground states of 32 S. Decays to excited states 
were identified by means of proton-7 coincidences with 
the 2230 keV and the 1548 keV 7 lines from the decay 
of excited states in 32 S. Figure [3] shows the 33 Ar /3-gated 
proton spectrum and the proton spectra coincident with 
the 7 lines at 2230 keV and the 1548 keV, respectively. 

At several instances, peaks due to proton emission to 
the ground state and to the excited states overlay. To 
ensure correct identification due to the different decay 
branches, the 7-gated proton spectrum for the 2230 keV 
7 line was normalized to the /3-gated spectrum by means 
of the proton peak at 1317 keV which is assumed to be 
only due to emission to the first excited state @. As 
this line is seen in the /3-gated spectrum and, of course, 
in the spectrum gated with the 2230 keV 7 line, both 
spectra can be normalized with respect to each other and 
a common relative normalization can be established. A 
fit of the 1317 keV peak in both spectra, 7 gated and 
singles, yielded a factor of 56 with which the 7-gated 
spectrum has to be multiplied to match both spectra. As 
the error of this number does not influence the final result 
in any significant way, it was not included. Figure[3J: and 
the shaded part of figure@]show this spectrum dominated 
by this proton peak. 

For the proton spectrum gated by the 1548 keV 7 line, 
we normalized the proton spectrum with the same factor 
as the spectrum in coincidence with the 2230 keV line and 
the ratio of the 7 detection efficiencies of the 1548 keV 
line and the 2230 keV line. This exploits the idea that in 
such an analysis the probability to observe protons in co- 
incidence with the 1548 keV line is only modulated with 
respect to the spectrum gated by the 2230 keV line by the 
7-detection efficiency The ratio of the two 7-ray efficien- 
cies was ei54s/e223o = 2.44%/1.89%. Thus, the spectrum 
generated with a 7 coincidence of the 1548 keV line was 
multiplied by this factor and the factor determined be- 
tween the /3-gated spectrum and the spectrum produced 
with the 2230 keV coincidence in order to match the rel- 
ative proton intensities for the ground state transitions, 
the transitions to the 2230 keV state and to the 3778 keV 
state in 32 S. 
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FIG. 3: (Color online) (a,b) The figure shows the /3-gated 
proton spectrum as obtained from the strips on the front side 
of detectors 1, 3, 5, and 6. The peak numbers correspond to 
those of table [TTJ Numbers below the curve are for protons 
decaying to the first excited state, but visible in the uncon- 
ditioned spectrum. The decay from the IAS in 33 CI to the 
ground state of 32 S is by far the most prominent peak, (c) 
The spectrum is generated from the same detectors, how- 
ever, with a condition that the 2230 keV 7 ray is observed 
in the germanium detectors, (d) The spectrum is obtained 
in a similar manner with a condition of the 1548 keV 7 ray. 
No subtraction is operated in the sense that spectrum (a,b) 
contains also proton groups decaying to the first excited state 
(e.g. the 1317 keV protons), whereas spectrum (c) contains 
also events in coincidence with the 1548 keV 7 line, as this 
2230 keV 7 ray always follows the 1548 keV 7 ray. The area 
of the proton peak at 3173 keV observed in the coincidence 
spectra with the 2230 and 1548 keV 7 lines gives the level of 
random coincidences of about 2xl0 -4 . For spectra (c) and 
(d) we indicate the reminder of this 3173 keV peak number 
21. 

To determine proton and 7-ray intensities, we used 
Gaussians for the peaks and straight-line fits for the 
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FIG. 4: (Color online) The figure shows the /3-gated proton 
spectrum as obtained from the strips on the front side of de- 
tectors 1, 3, 5, and 6 overlaid with the spectrum obtained 
from the 2230 keV 7-ray coincidence spectrum. The normal- 
ization of the second spectrum is such that the intensities of 
the 1317 keV peak, being to 100% due to the 2230 keV 7 co- 
incidence, match. To generate the 7-coincident spectrum, a 
spectrum was first generated with a 20 keV wide condition on 
the 2230 keV 7 line. From this spectrum, we subtract then a 
similar spectrum conditioned by a gate left and right of the 7 
line. 



background. All intensities were first normalized with 
respect to the most intense peak (the 3173 keV proton 
line and the 810 keV 7 line, respectively) to obtain rela- 
tive branching ratios. 

The relative intensities of the proton groups can eas- 
ily be converted into absolute intensities by assuming 
that the absolute intensity of the feeding of the IAS at 
5549 keV can be reliably determined. For this purpose, 
one usually assumes that the Fermi strength can be cal- 
culated more or less model independently b y n eglecting 
Coulomb and radiative corrections (see e.g. jT^| ) . In the 
case of 33 Ar as for other odd-mass nuclei, the additional 
problem of a GT contribution to the decay to the IAS 
arises. This contribution can either be determined ex- 
perimentally (see e.g. [§, HU), be neglected assuming its 
contribution to be small, or taken from theoretical pre- 
dictions. We will calculate the Fermi part neglecting the 
corrections mentioned above and take the GT contribu- 
tion from Ref. 



11] 



The Fermi strength can be calculated from the relation 
between the ft value and the Fermi matrix element Mp: 
ft = Kj{g' 2 v * Mp), where K is a constant and g\ is the 
effective vector coupling constant of the weak interaction. 
For K/g^, we use a value of 6144.2(13) s deduced from 
the super-allowed + — > + decays [lg. The matrix ele- 
ment for nuclei with an isospin quantum number T = 3/2 
is y/3. This allows the Fermi ft value to be calculated to 
be 2048 s. 
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Peak 


Energy 


BR (%) 


BR (%) 


BR (%) 


BR (%) 


number 


(kcV) 


this work 


Schardt ct al. [91 


Borgc ct al. [8 


average 


protons to the ground state of 32 


S 






3 


1645(2) 


0.411(20) 


0.391(6) 


0.343(10) 


0.380(16) 


7 


1781(2) 


0.471(22) 


0.459(6) 


0.434(10) 


0.453(8) 


9 


2100(3) 


2.73(12) 


2.37(2)+0.35(l) W 


2.37(2) 


2.54(12) 


12 


2481(2) 


0.362(17) 


0.353(6) 


0.333(10) 


0.349(7) 


14 


2744(3) 


0.0483(44) 


0.0403(31) 


0.0454(50) 


0.0435(25) 


17 


2941(4) 


0.0748(55) 


0.0713(31) 


0.1222(81) 


0.0772(107) 


21 


3173(3) 


31.0(14) 


31.0(14) 


31.0(14) 


31.0(14) 


22 


3350(4) 


0.0918(48) 


0.0310(62) 


0.7573(404) 


0.0753(496) 


26 


3857(3) 


0.735(34) 


0.716(6) 


0.808(20) 


0.724(18) 


32 


4719(5) 


0.00079(10) 


- 


- 


0.00079(10) 


33 


4860(4) 


0.0097(8) 


0.0152(31) 


0.0232(30) 


0.0108(24) 


35 


5039(4) 


0.224(12) 


0.217(6) 


0.333(10) 


0.245(34) 


36 


5101(4) 


0.0470(50) 


0.0589(31) 


0.1111(101) 


0.0592(103) 


37 


5225(4) 


0.0234(24) 


0.0288(16) 


0.0545(40) 


0.0298(59) 


38 


5317(3) 


0.0083(12) 


0.0133(12) 


0.0182(30) 


0.0113(22) 


39 


5623(3) 


0.124(7) 


0.118(25) 


0.222(10) 


0.155(32) 


10 


5723(3) 


0.092(5) 


0.062(15) 


0.172(10) 


0.103(22) 


11 


5855(9) 


0.00284(40) 


- 


0.01212(202) 


0.0032(18) 


42 


6011(10) 


0.00100(15) 


- 


- 


0.00100(15) 


43 


6100(10) 


0.0138(15) 


- 


0.0212(40) 


0.0147(24) 


41 


6344(8) 


0.00053(9) 


- 


0.00505(202) 


0.00054(21) 


45 


6389(10) 


0.00027(8) 


- 


- 


0.00027(8) 


46 


6480(10) 


0.0103(11) 


- 


0.0162(20) 


0.0116(24) 


47 


6628(10) 


0.00170(23) 


- 


0.00404(202) 


0.00173(26) 


48 


6657(9) 


0.00049(10) 


- 


- 


0.00049(10) 


49 


6715(9) 


0.00100(12) 


- 


0.00364(162) 


0.00102(20) 




6800 - 6900 


0.00023(9) 


- 


- 


0.00023(9) 




6900 - 7000 


0.00007(3) 


- 


- 


0.00007(3) 




7000 - 7100 


0.00032(4) 


- 


- 


0.00032(4) 




7100 - 7200 


0.00000(0) 


- 


- 


0.00000(0) 




7200 - 7300 


0.00012(3) 


- 


- 


0.00012(3) 




7300 - 7400 


0.00010(3) 


- 


- 


0.00010(3) 




7400 - 7500 


0.00008(3) 


- 


- 


0.00008(3) 




7500 - 8000 


0.00006(3) 


- 


- 


0.00006(3) 




8000 - 9000 


0.00004(3) 


- 


- 


0.00004(3) 


protons to the first excited state 


of 32 S 






1 


762(10) 


0.0202(17) 


- 


- 


0.0202(17) 


2 


1317(8) 


0.168(9) 


0.180(3) 


0.191(8) 


0.180(4) 


3 


1691(6) 


0.0332(32) 


0.0319(16) 


0.0464(61) 


0.0329(22) 


6 


1764(5) 


0.0081(13) 


0.0220(31) 


- 


0.0102(49) 


8 


2024(5) 


0.0043(7) 


- 


- 


0.00428(71) 


11 


2370(5) 


0.0153(12) 


0.0158(31) 


0.0192(30) 


0.0159(9) 


13 


2710(7) 


0.0069(12) 


- 


- 


0.0069(12) 


15 


2810(10) 


0.00141(14) 


- 


- 


0.00141(14) 


16 


2886(7) 


0.0376(35) 


0.0341(31) 


0.0646(61) 


0.0394(70) 


18 


2957(7) 


0.0359(32) 


0.0434(31) 


- 


0.0397(38) 


20 


3066(6) 


0.00440(61) 


- 


0.7068(2020) 


0.0045(20) 


23 


3469(6) 


0.0531(40) 


- 


0.4140(404) 


0.057(36) 


24 


3515(6) 


0.0150(25) 


0.0065(16) 


- 


0.0089(38) 


25 


3576(5) 


0.0085(16) 


- 


0.1010(101) 


0.011(14) 


27 


3926(5) 


0.0082(13) 


- 


0.0192(40) 


0.0091(31) 


28 


4209(5) 


0.00645(83) 




0.0172(30) 


0.0072(27) 


29 


4330(8) 


0.00142(40) 




0.0081(20) 


0.0017(13) 


30 


4474(5) 


0.00367(55) 




- 


0.00367(55) 


31 


4505(6) 


0.00467(62) 




0.01820(51) 


0.0049(16) 


34 


4923(6) 


0.00066(16) 






0.00066(16) 


protons to the second excited state of 32 S 






4 


1665(6) 


0.0060(11) 


0.0099(16) 




0.0074(19) 


10 


2368(6) 


0.0012(3) 






0.0012(3) 


19 


3016(10) 


0.0007(2) 






0.0007(2) 



TABLE II: Summary of proton energies in the laboratory frame (from this work) and absolute branching ratios from the 
present work, from Schardt et al. [§] and from Borge et al. §. Eight proton lines attributed by Borge et al. to a transition 
to the ground state were newly attributed to proton emission to the first excited state. One proton group of Schardt et al. is 
attributed to the decay to the second excited state. Two proton lines observed by Schardt et al. [3] are neither resolved in 
the present work nor in the work of Borge et al. [8|. The absolute intensities have been determined by normalizing all proton 
group with respect to the 3173 keV line for which the absolute branching ratio was established as described in the text. The 
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We take the GT to Fermi ratio from Ref. where 
it was determined to be 0.044±0.002. Adding this value 
to the Fermi matrix element and taking into account the 
ratio between the squares of the vector and the axial- 
vector coupling constants of 1.2695(29) [l7[, we derive 
a total ft value for the decay to the IAS of 1962(5) s. 
As we have neglected the corrections mentioned above, 
we add an additional uncertainty of 1% for this number 
which yields the final ft value for the IAS of 1962(20) s. 

this work Borge et al. f8] average 

E-y (keV) BR (%) (keV) BR (%) BR (%) 



33 C1 810.6(2) 
33 C1 1541.4(6) 
32 S 2230.3(5) 
33 C1 2352.5(6) 
33 S 2866.2(9) 



100(1) 
3.6(2) 
6.1(2)* 
1.3(2) 
0-6(1) 



810.3(5) 
1541.5(5) 
2230.6(9) 
2352.2(9) 
2867.7(9) 



100(3) 
2.4(5) 
1.7(5) 
1.7(5) 
1-3(1) 



100 
3.4(4) 

1.4(2) 



TABLE III: Gamma-ray energies and their relative branching 
ratios from the present work and the work of Borge et al. [8J . 
Averages are given only for 7 rays following the decay of 33 Ar. 
For the 7 ray at 2230 keV, the relative branching ratio in our 
experiment indicated by the star is too high, when it is cal- 
culated with respect to the 810 keV line. The former 7 ray is 
always preceded by a proton which triggers the DSSSDs, while 
for the latter the P particles have to trigger. Evidently, the 
P trigger probability is much smaller than the proton trigger 
efficiency which then yields higher relative branching ratios 
for /3p7 events. The literature value for the last 7-ray energy 
is 2867.72(2) keV. 

To determine the absolute branching ratio to the 
IAS, we now need the statistical rate function or 
Fermi function calculated from the /3-decay Qec value 
(11.6193(6) MeV [llj]) from which we subtract the excita- 
tion energy of the IAS (5.5478(9) MeV) Q3] and the half- 
life of 33 Ar. For the half-life, we use the average of differ- 
ent literature values [U-IH for which we find 173.9(9) ms. 
The statistical rate function is calculated to be 3501(4). 
With these numbers, we determine a branching ratio to 
the IAS of 31.0(14)%. 

The branching ratio to the ground state of 33 CI could 
not be determined in our experiment. We therefore as- 
sume that mirror symmetry is a good approximation and 
calculate its feeding from the ft value of the mirror de- 
cay. 33 P decays with a 100% branch to the ground state 
of 33 S. The half-life (25.34(12) d) and the Q E C value 
(248.5(11) keV) [l8[ enables us to determine the ft value 
to be 107280(2250) s (logft = 5.03(1)). Using this ft 
value for the decay to the ground state of 33 CI, we deter- 
mine a branching ratio of 18.7(4)%. 

The difficulty is now to convert the branching ratios for 
the feeding of bound levels in 33 CI to absolute branch- 
ing ratios. We will test two ways of determining these 
branching ratios. In the first procedure, we use the to- 
tal proton branching ratio and the branching ratio to the 
ground state to deduce the feeding of bound excited lev- 
els as the three contributions should add up to 100%, as 
long as 7 decay of proton-unbound levels is negligible. 




750 
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FIG. 5: (Color online) Gamma-ray spectrum registered with 
the three germanium clover detectors. The peaks labeled 
belong directly or indirectly to the decay of 33 Ar. The in- 
set shows the region around the 1541 keV and 1548 keV 7 
rays. The 1548 keV peak from the decay of the second ex- 
cited state in 32 S is very weak and only visible as a shoulder 
of the stronger 1541 keV peak from the decay of the sec- 
ond excited state of 33 CI. The relative intensities of all rel- 
evant peaks can be found in table [HI] Background lines are 

from 60 Co (1173 keV, 1332 keV), from a contamination of the 
beam with 35 Ar (1220 keV, 1764 keV, 2694 keV), from 40 K 
(1461 keV), and from the decay of the first excited state of 
208 Pb (2615 keV). 



In the second, we use the decay of 33 CI and the relative 
intensities of the 7 transition in 33 CI to determine the 
feeding of bound levels. 

For the first procedure, we can obtain the total pro- 
ton branching ratio in two ways: i) by summing all pro- 
ton branching ratios in table HH ii) by integrating the 
/3-gated proton spectrum assuming that it contains all 
proton events. However, in second this case, one has to 
cut the (3 particles by means of a threshold which we 
chose, somewhat arbitrarily, to be 1000 keV. We note 
that the final result changes by 0.8%, if we integrate the 
proton spectra starting at 900 keV. We include this in 
the final error. The total proton rate can then be nor- 
malized by means of the known absolute branching ratio 
of the IAS of 31.0(14)%. We determine a ratio between 
the number of counts for the IAS proton line and the to- 
tal proton intensity of 78(2)%. In the former procedure, 
where we sum the branching ratios from table HJ we ob- 
tain a total proton branching ratio of 36.9(17)%. The 
latter procedure gives a value of 39.7(12)%. We adopt a 
final value of 38.8(13)% for this first procedure. 

With the branching ratio to the ground state of 33 CI 
of 18.7(4)% and the total proton branching ratio of 
38.8(13)%, we derive a branching ratio to bound excited 
states in 33 C1 of 42.5(17)%. From this number and the 
relative 7 branching ratios (see figure [5] and table IIII|) , 
one can determine the absolute branching ratios for the 
feeding of the 810 keV level of 40.5(16) % and for the 
2352 keV level of 2.0(3)%. 
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TABLE IV: Excitation energies, 7 and laboratory proton energies, branching ratios (averages from the present work, Borge et 
al. 0, and Schardt et al. 0), individual B(GT) values and the summed B(GT). In cases, where the excitation energy is missing 
for the present work, the proton is emitted from the same level as the proton in the previous line. Proton energies labeled with 
# are due to protons populating the first excited state, whereas proton energies labeled with are due to protons populating 
the second excited state. Shell-model values for the excitation energies and the branching ratios are given for levels below the 
IAS. They were obtained using the USD effective interaction [2l|. The excitation energies determined in the present work are 
„i„„ ,,„!,,„„ T?„,i+ ForSl ;r +u„ !„,,„! „ T „„ j„+„^™;„„j f„„„, a J„„„,, \\r„ +u„ „„™„ 
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In the second procedure to determine the total branch- 
ing ratio for decays to bound states in 33 CI, we calculate 
these branching ratios from the observation of the decay 
of 33 C1. The 7 ray at 2868 ke V from this decay has a 
branching ratio of 0.44(6)% 22]. If we use the observed 
relative branching ratio for this 7 ray from table Mil and 
the possible branches populating the ground state of 33 CI, 
we determine a branching ratio of the 810 keV 7 ray of 
33 C1 of 51.4(120)% and of 2.6(6)% for the branching ra- 
tio to the 2352 keV level. These values agree with those 
determined with the first procedure, but they are signif- 
icantly less precise. We therefore adopt the values ob- 
tained with the first procedure. 

The procedures just laid out assume that there is no 
7 decay of proton decaying levels. In addition, forbid- 
den transitions could falsify the GT strength distribu- 
tion calculated in the following. The possible presence of 
forbidden transitions is shortly discussed below. 

With these absolute branching ratios for the feeding 
of the ground state of 33 CI as well as of the bound and 
unbound excited states of 33 C1, one can determine the 
feeding of all states in the decay of 33 Ar, the Gamow- 
Teller strength B(GT) and its sum. These data are given 
in Table [TV] For the 762 keV and 2941 keV proton lines, 
we adopted the level energy calculated only from the 
2941 keV protons, as the low-energy protons are much 
stronger affected by the energy loss in the silicon dead 
layer and in the mylar foil of the catcher than the high- 
energy protons. For all other decays, we use average val- 
ues from decays to the ground and first (second) excited 
states. In addition, we determine the log(ft) values for 
all observed transitions which are given in figure [9] 



IV. DISCUSSION OF THE RESULTS 
A. Comparison with previous experimental work 

In the work of Borge et al. H , the energy calibration 
of the silicon detector at higher energies suffers from the 
fact that the calibration was performed with a particles 
and that the pulse height defect [23[ was not corrected 
for. Therefore, with increasing energy above the IAS, the 
proton energies determined become smaller and smaller. 
In order to find the proton lines in Borge et al. corre- 
sponding to proton lines of the present work, we corrected 
for this effect of the Borge et al. data. 

Compared to previous work (8Hlfj|. we observe all the 
proton lines identified in the work of Borge et al. Q. 
However, eight of them are now attributed to proton 
emission to the first excited state rather than to the 
ground state (see table |n|. Generally, the branching ra- 
tios of Borge et al. are higher than our values. As for 
Schardt et al. |9|, the proton group at 1750 keV corre- 
sponds most likely to our proton energy of 1764(5) keV 
observed in coincidence with the 2230 keV 7 ray. The 
lines at 2121 keV and 2096 keV observed by Schardt et 
al. are not resolved in the present work. They corre- 



spond to our line at 2100(3) keV. The 1665 keV line is 
assigned to the decay to the second excited state, as ob- 
served in the present work. Honkanen et al. fiol ] identified 
two weak lines at 3485 keV and 5658 keV. The first one 
corresponds probably to our line at 3469(6) keV which 
we observe in coincidence with the 7 ray at 2230 keV (see 
below). For the higher-energy line, we can not exclude 
that it corresponds to our proton line at 5623(3) keV. 

Table ITll gives the proton group energies, their absolute 
intensities from the different experiments and the average 
intensities. The branching ratios of the previous work 
were updated with respect to the branching ratio to the 
IAS as determined above. For the experimental proton 
energies, we give only our energies. 

In table IIV1 we compare the levels observed in the 
present work with those found in the latest review of 
properties of the mass A=33 chain by Endt [20]. The 
levels given in this compilation come both from /3-decay 
studies [§[ and from reaction work [24Tj27| . For levels de- 
termined by different experimental approaches, we find 
reasonable agreement. For higher-energy levels, we just 
compare with those deduced from reaction work as the 
comparison with /3-decay work has already been done in 
table HIl Therefore, in table HVl we only compare with 
previously established states deduced from work other 
than /3 decay. Similarly, the spin values shown come only 
from reaction studies. 

B. Comparison with shell- model calculations and 
the GT strength distribution 

If we compare our results to shell-model calcula- 
tions [l!| (see table IIV[ ). we find that below the IAS only 
one predicted state at 3.863 MeV with a branching ratio 
of 0.0747% is not observed experimentally. All other pre- 
dicted states are observed experimentally with branching 
ratios in most cases close to the predicted ones. 

For a comparison of the experimental GT strength 
distribution with theoretical predictions, the theoretical 
strength has to be quenched. The origin of this quenching 
is not completely clear, but excitations of sub-nucleonic 
degrees of freedom, notably of the A resonance, and, 
probably to a much larger degree, contributions from 
outside the considered model space are most often sus- 
pected (see e.g. [H). A typical quenching factor for the 
sd shell is q 2 = 0.5. With our data, we can check this 
quenching factor. We will use experimental B(GT) val- 
ues deduced for individual transitions below the IAS and 
compare them to the B(GT) strength for the equivalent 
transitions (see table IIV[) as predicted by theory. This 
comparison is shown in figure [51 A fit of these theoreti- 
cal data as a function of the experimental data yields a 
quenching factor of 0.49(4). This is in perfect agreement 
with the generally accepted quenching in the sd shell. 

Using this quenching factor, we can compare the sum- 
med experimental B(GT) strength with the prediction 
from the shell model. For this comparison, we use cal- 
culations with three different effective shell-model in- 
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FIG. 6: (Color online) Comparison of the GT strength be- 
tween experiment and shell-model calculations for eight in- 
dividual transitions of 33 Ar for which a correspondence be- 
tween the shell-model calculations with the USD interaction 
and experiment was established in table HVl The fit yields a 
quenching factor of 0.49(4). 

teractions for the sd model space: the USD interac- 
tion [2l| and two newly determined interactions USDa 
and USDb [28|. For this purpose, we convert the exper- 
imental branching ratio into a B(GT) strength for each 
individual level observed. As a function of the excita- 
tion energy, this strength is then summed. The result 
is shown in figure [7] At low excitation energies, almost 
perfect agreement of the experimental data with all three 
interactions is obtained. At higher excitation energy, the 
new interactions USDa and USDb seem to slightly better 
fit the data. Nonetheless, for all three interactions, the 
agreement is remarkable. 

The procedure just described has one drawback: Only 
decay strength identified as a peak usually is taken into 
account. Very weak decay branches are omitted. There- 
fore, we developed a procedure which uses all strength 
present in the proton spectra. For this purpose, we take 
the total charged-particle spectrum as obtained by means 
of the silicon detectors and subtract first the /3 back- 
ground by a fit of the low-energy part of this spectrum. 
In a second step, we subtract the decay strength to the 
first and second excited state. These contributions are 
obtained from the proton spectrum conditionned by the 
detection of these 7 rays renormalized by the gamma ef- 
ficiency. This yields the proton decay strength to the 
ground state. Now we have to convert this spectrum bin 
by bin to a center-of-mass (CM) spectrum. To this we 
add the decay spectrum (again converted to the center- 
of-mass) to the first excited state shifted by the excita- 
tion energy of this state. Similarly, we treat the spec- 
trum for decays to the second excited state. After a 
correction for the proton separation energy, this yields 
the excitation energy spectrum to which we have added 
the decay to bound states (see figure [5^) . After calcu- 
lating the statistical rate function / for each energy bin, 
we can determine the B(GT) value for each energy bin 
(figure |5}d). This spectrum can be compared to the theo- 
retical B(GT) values for individual transitions (figure [8fc) 
and the summed B(GT) distribution (figure [SJi). At high 
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FIG. 7: (Color online) Summed Gamow- Teller strength distri- 
bution for the decay of 33 Ar. The experimental distribution is 
compared to shell-model calculation by B.A. Brown [19] with 
the universal sd shell interaction USD and two updated ver- 
sions of the same interaction [28] . For these calculations, we 
used a quenching factor of 0.5, the generally accepted value. 
The experimental data correspond only to resol ved proton 
groups obtained from the average values of Table IIVI . 

excitation energies, the summed experimental B(GT) dis- 
tribution exceeds the theoretical one with a quenching of 
0.5. If we integrate the B(GT) up to 11 MeV we get 5.3 
which correspond to a quenching of 0.59. To discuss the 
contribution at very high excitation energy is always del- 
icate, enough to say that if we integrate in the full QeC 
window of 11.6193 MeV the cumulated B(GT) grows to 
41.2. 

A possible explanation for this large B(GT) sum is 
that we shift too much strength to high excitation ener- 
gies by attributing too much decay strength to decays to 
excited states, in particular to the second excited state. 
To test this assumption, we calculated the GT strength 
distribution also by omitting the second excited state and 
both excited states. However, even when omitting both 
excited states, we still exceed the quenched theoretical 
B(GT) values. 

Another possible explanation could be that we deal 
to some extent with forbidden transitions. In the latest 
analysis of Singh et al. [29j, the log(ft) values of allowed 
transitions are distributed over a large range from 3 to 7, 
whereas forbidden transitions start with values of about 
5. As can be seen in figure HO we determined indeed 
high log(ft) values which could be due to forbidden tran- 
sitions. However, as we have no means to distinguish 
between allowed and forbidden transitions, we assume 
that all transitions are allowed, but we give their spins 
in parenthesis to indicate this uncertainty. 

Finally, background is of course also a concern. 
Proton-proton pile up might create high-energy events 
which get an important weight due to the larger space 
phase factor. However, we do not have any means to 
distinguish such possible background counts from real 
counts. 

The good agreement between experimental and theo- 
retical B(GT) strength up to rather high excitation ener- 
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FIG. 8: (Color online) (a) Proton and 7-ray spectrum trans- 
formed into a spectrum of excitation energies. The three tran- 
sitions to bound states are the discrete lines at low energies. 
The solid line shows the statistical rate function (arb. units), 
(b) The spectrum gives the 7 and proton transitions converted 
into the Gamow- Teller strength of the individual experimental 
transitions (see text for details) . The contribution of the IAS 
is removed, (c) Gamow- Teller strength for individual transi- 
tions as calculated by means of shell-model calculations with 
the USD interaction, (d) Integrated Gamow- Teller strength 
as a function of the excitation energy as determined from the 
proton spectrum by adding the three transitions to bound 
states for the experimental spectrum (full line) and the shell- 
model calculations with the USD interaction (dashed line). 



gies shows that the quenching is the same at low and at 
high excitation energies. A comparison of figures [7] and [5] 
shows that no major proton groups are missing in ta- 
ble |TT] and basically all strength is identified in the peaks. 
If this were not the case, the second procedure arriving 
at figure |8] should yield more B(GT) strength than the 
procedure which produces the spectrum of figure [7J 



C. Spins and parities 

Under the assumption of allowed transitions, all lev- 
els populated in the /3 decay of 33 Ar (r T = l/2 + ) have 
spin/parity l/2 + or 3/2 + . However, certain characteris- 
tics observed in the present study may allow to distin- 
guish between these two possibilities. We will use the 
observed proton lines to propose spins for some levels. 
This will be discussed in the following. 

Proton emission from a nuclear state is governed by 
two effects: the barrier penetration and the spatial over- 
lap of the initial and the final states. For a given nucleus, 
the barrier penetration depends on the energy available 
for the proton and the barrier height. The contribution 
of the Coulomb barrier is the same for all states and 
all emissions. However, for a transition with a non-zero 
angular momentum, the angular momentum barrier is 
added. Therefore, proton transitions to the ground state 
are favored by the available energy, except if angular mo- 
mentum has to be " carried away" by the proton. If states 
emit protons only to an excited state in the proton daugh- 
ter nucleus, this can be an indication that this transition 
has a lower angular momentum than the ground state 
transition, a fact which would enable us to distinguish 
between two possible spin assignments. 

Such an assignment, however, neglects the possibility 
that strongly different overlaps of the initial and final 
wave functions can have a similar effect. In principle, 
model calculations would have to be performed for all 
states involved before such an assignment can be made. 
We believe that this is not very meaningful, as the states 
calculated e.g. in shell-model calculations might differ 
strongly from the observed states, even if the excitation 
energy is equivalent. 

We assign tentatively a spin of 3/2 + to all states 
which decay only to the first excited state in the proton- 
daughter nucleus 32 S. We give these spin/parities in 
parenthesis to make clear that this is a tentative assign- 
ment. 

The states at 7762 keV and at 8491 keV decay to the 0+ 
ground state and the + second excited state. This indi- 
cates that the emitting state has most likely a spin/parity 
of 1/2+ which yields an I — proton emission. 

For other states, barrier penetration calculations (see 
tabic |V} may elucidate the situation. The barrier pen- 
etration is calculated from Coulomb wave functions [30l ] 
which yields partial half-lives for both spin possibilities 
for decays to the ground and excited states. The ratios 
of the partial half-lives for emission to the ground and 
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excited states for both spins are then compared to the 
ratio of the experimental branching ratios. If we believe 
that the situation is relatively clear and one spin assign- 
ment gives a much closer agreement than the other spin, 
we propose a spin. The procedure used is tested with 
and works well for the state at 7292 keV, for which we 
find a spin of 3/2+ in agreement with the literature [20| . 

For the state at 9152 keV which decays to the three 
lowest states of 32 S, we determine that for a 3/2 + emit- 
ting state the decay to the first excited state is indeed 
the fastest transition thus receiving the highest branch- 
ing ratio. For a spin 1/2+, the decay to the ground state 
is calculated to be a factor of ten faster than the de- 
cay to the other two states, opposite to the experimental 
observation. We therefore attribute a spin 3/2 + to this 
state. 

We underline again that the present spin assignments 
have to be confirmed by other means before they should 
be accepted. Therefore, they appear in parenthesis in 
figure [9] 



D. Search for candidate states for isospin mixing 

Proton emission from the T=3/2 IAS in 33 CI can only 
take place due to isospin impurity either in the proton 
emitting state or in the proton daughter state. Therefore, 
Honkanen et al. [l(| searched for states in the vicinity of 
the IAS with which the IAS can mix. We identified three 
of the four states discussed by Honkanen et al. There is 
one state about 230 keV below the IAS and two other 
states 180 keV and 310 keV above the IAS. However, we 
do not have evidence for a state about 100 keV below 
the IAS which would be the best candidate for mixing 
between such a state (T=l/2, F=l/2+) and the IAS 
(T=3/2, F=l/2+). 

The proton group identified by Honkanen et al. [l(| at 
3073 keV corresponds to our peak at 3066 keV. However, 
we clearly identify the proton group to decay to the first 
excited state and not to the ground state as suggested 
by Honkanen et al. Therefore, the emitting state is no 
longer in the vicinity of the IAS. 

We do not completely exclude the possibility of such 
a close-by state, as our resolution is maybe not good 
enough to resolve several small peaks in the close vicinity 
of the strong IAS. However, if such a state is absent, 
the mixing has to take place with states further away 
from the IAS. A better resolution experiment with higher 
statistics is certainly needed to answer this question. 

A state with (T=l/2, F=l/2 + ) about 100 keV be- 
low the IAS was also observed in proton scattering ex- 
periments on 32 S [25j. This would of course be a good 
candidate for isospin mixing with the IAS. As mentioned 
above, we do not observe such a state. However, we can- 
not exclude that our proton peak at 3066 keV contains 
also a small branch to the ground state of 32 S. 

It is interesting to note that the shell model using the 
three different USD interactions @, [28| does not predict 



any (T=l/2, F=l/2+) state close to the IAS. The states 
close to the IAS have all a spin of I 7r =3/2 + . The clos- 
est I 7r =l/2+ state is at about 400 keV above the IAS for 
the USD interaction, 550 keV above for the USDa in- 
teraction, and 450 keV above the IAS for USDb. Below 
the IAS, the first (T=l/2, F=l/2+) state is typically 
800-900 keV away from the IAS. 

E. Decay scheme of 33 Ar 

The information accumulated in the present paper to- 
gether with the results from previous work enables us to 
establish a rather complete decay scheme of 33 Ar with 
decays to bound states and proton unbound states which 
then decay to the three lowest states of 32 S by proton 
emission. Figure [5] gives this decay scheme. We remind 
the reader that the branching ratios for the three proton- 
bound states are based on the assumption of mirror sym- 
metry of the ground states of 33 Ar and 33 P. 

V. SUMMARY 

Using a novel experimental setup of high granularity 
and high efficiency for charged-particle detection com- 
bined with a high-efficiency 7 array allowed a complete 
study of the decay of 33 Ar to be performed. The experi- 
mental results are to a large extent in agreement with pre- 
vious experimental studies of this nucleus. They enabled 
us to establish a quasi complete decay scheme of 33 Ar and 
to compare the experimental results with shell-model cal- 
culations using different effective interactions optimized 
for the sd shell-model space. The excellent agreement 
obtained testifies of the high quality of the shell model 
in this region of the chart of nuclei. The comparison 
of experimental B(GT) strength to well identified levels 
and theoretical calculations yields a quenching factor of 
0.49(4), in excellent agreement with the accepted value. 
States in close vicinity of the isobaric analog state of the 
33 Ar ground state in 33 CI are searched for and discussed 
in terms of isospin mixing. However, a state identified 
before could not be observed in the present work. 
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TABLE V: Barrier penetration calculations used to tentatively assign spins to certain levels. The first column gives the 
excitation energy of the state, followed by the two or three proton energies from this state. The next two columns give the 
angular momentum used in the calculations which is deduce from the two spin possibilities of the proton emitting state, 1/2 or 
3/2, given in the following column. The next column is the barrier penetration half-life calculated with the angular momentum 
1=0 or 1=2. The sixth column is the experimental branching ratio of the two/three proton transitions from the state. These 
branching ratios are given only once for each proton group. For some states, in the final column, we give our tentative spin 
assignments. 
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FIG. 9: (Color online) Decay scheme of 33 Ar. All the decay branches, either to particle-bound states or to proton-unbound 
states decaying by proton emission to the three lowest states of 32 S, are shown. The excitation energy of the decaying state, 
the branching ratio and the log(ft) value are given. Spins tentatively assigned in the present work to be l/2 + or 3/2 + are given 
in parentheses. The spins of levels for which both spins (l/2 + , 3/2 + ) are given in parenthesis are determined from /3 decay 
selection rules assuming allowed transitions. Other spin/parities are from Ref. [28] . 
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